A generalized thermodynamics for power-law statistics 
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We show that there exists a natural way to define a condition of generalized thermal equilibrium 
between systems governed by Tsallis thermostatistics, under the hypotheses that i) the coupling 
between the systems is weak, ii) the structure functions of the systems have a power-law dependence 
on the energy. It is found that the q values of two such systems at equilibrium must satisfy a 
relationship involving the respective numbers of degrees of freedom. The physical properties of a 
Tsallis distribution can be conveniently characterized by a new parameter r\ which can vary between 
and +00, these limits corresponding respectively to the two opposite situations of a microcanonical 
distribution and of a distribution with a predominant power-tail at high energies. We prove that the 
statistical expression of the thermodynamic functions is univocally determined by the requirements 
that a) systems at thermal equilibrium have the same temperature, b) the definitions of temperature 
and entropy are consistent with the second law of thermodynamics. We find that, for systems 
satisfying the hypotheses i) and ii) specified above, the thermodynamic entropy is given by Renyi 
entropy. 

PACS numbers: 05.70.-a, 05.20.-y, 05.90.+m 
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I. INTRODUCTION 



In recent years a considerable amount of theoretical 
research activity [E El H, EL IE @ nas been stimulated by 
the conjecture that a certain class of dynamical systems 
might be governed by a peculiar thermal statistics, in 
which the usual Boltzmann-Gibbs exponential distribu- 
tion is replaced by a type ofpower-law distribution orig- 
inally proposed by Tsallis [7[. Such a conjecture seems 
to be supported by phenomenological and computational 
evidence coming from various domains of physics and 
other disciplines. 

It has been proposed in various ways to define a gen- 
eralized entropy function, according to formulae involv- 
ing a parameter q directly related to the characteristic 
exponent of the power-law distribution (whence the fre- 
quently used name of "^-distributions" ) . It appears then 
natural to extend to these systems also the concept of 
temperature, in such a way that the validity of the sec- 
ond law of thermodynamics is preserved. However some 
difficulties have to be faced in order to define the condi- 
tion of thermal equilibrium and to consistently general- 
ize for these systems the zeroth law of thermodynamics. 
The simplest possibility is apparently that of treating 
systems at thermal equilibrium as statistically indepen- 
dent of one another [8j|. However, it has already been 
recognized that the assumption of such an independence 
has no fundamental justification. Furthermore, the prob- 
ability distribution for a composed system obtained by 
multiplication of two power-law distributions does not 
belong in general to the same category. This problem is 
not present in the Boltzmann-Gibbs case thanks to the 
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peculiar property of the exponential function. As a way 
out of this difficulty, a generalization of the zeroth law in 
the spirit of superstatistics [9( has recently been proposed 

A different approach to the thermodynamics of sys- 
tems with power-law statistics has been developed fol- 
lowing the observation that a distribution of this type 
can also be formally derived using standard arguments 
of equilibrium thermodynamics. In fact, in the same 
way as the usual Boltzmann-Gibbs distribution is de- 
rived for a system in thermal contact with an infinite 
heat bath, a power-law distribution is obtained when 
the heat bath is instead characterized by a finite number 
of degrees of freedom and a power-law density-of-states 
function [111 |12| . 1 1 31 ] . In such a situation, it becomes nat- 
ural to define the physical temperature by making use of 
the equipartition theorem [l4j]- It has been shown that, 
when the inverse of this temperature is used as the inte- 
grating factor of the exchanged heat, in accordance with 
the second law of thermodynamics, the definition of the 
entropy has to be modified with respect to the original 
Tsallis formula [HI, [l6| . It has however to be noted that 
the presence of a finite heat bath can only be assumed 
when the power-law distributions present a finite energy 
cut-off, which is not the case for all the possible ranges 
of the parameter q. 

The distinction between "observed system" on the one 
side, and "heat bath" on the other, which has been fol- 
lowed in Refs. 

mi mm mm 

in analogy with famil- 
iar considerations of standard thermodynamics, appears 
somewhat fictitious when both systems are supposed to 
be finite. In the present paper we shall then abandon 
the concept of heat bath altogether, and we shall only 
use that of mutual thermal (mcta)equilibrium, in a form 
which can in principle be applied to systems obeying any 
(a priori arbitrary) generalized type of thermal statistics. 
If we consider a system composed of parts in mutual ther- 
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mal contact, it is obvious that when the system is in a 
stationary state, its parts are in thermal metaequilibrium 
with one another (of course, the parts in question may 
be either of homogeneous or of heterogeneous physical 
nature). We shall investigate, with the widest possible 
generality, the logical consequences of this simple fact in 
connection with the laws of thermodynamics. In partic- 
ular, in order to determine the correct definition for the 
absolute temperature, we shall naturally impose the con- 
dition that, in accordance with the zeroth law, systems at 
equilibrium must have the same temperature. We shall 
see that such an approach applies equally well to Tsallis 
distributions having a finite energy cutoff, and to those 
showing instead a power tail extending to infinitely large 
energies. 

From a mathematical point of view, in order to apply 
in a consistent way the concept of thermal equilibrium to 
states described by Tsallis statistics, we have first of all 
to establish under which conditions the probability dis- 
tribution for a composed system, and the marginal distri- 
butions which are derived from it, can all have simulta- 
neously a power-law dependence on the energy. We will 
show that a sufficient condition is that also the densities- 
of-states of all the subsystems considered are given by 
power functions of the energy. Of course we know that 
this is indeed the case for ideal gases, and represents a 
good approximation for a quite large class of systems. For 
these systems we can then interpret Tsallis statistics as 
a general description of metaequilibrium states, which is 
consistent with the zeroth law of thermodynamics. With 
respect to standard statistical mechanics [l?], [3 , which 
is only based on the Boltzmann-Gibbs distribution, the 
new theory contains a set of additional parameters, the 
q values of the systems considered. It will be shown that 
these parameters, in order for thermal equilibrium to es- 
tablish, must be in a certain relationship with the respec- 
tive numbers of degrees of freedom. 

The first and second laws of thermodynamics maintain 
their validity for systems governed by power-law distri- 
butions. An important point, which will emerge from our 
analysis, is that these laws, when considered in combina- 
tion with the zeroth one, univocally determine the ex- 
pression of both the temperature and entropy functions. 
It will turn out that, for systems to which our starting 
hypotheses can be applied, the correct expression for the 
entropy coincides with that originally proposed by Renyi 

Mm 

The consideration of power-law densities of states will 
also allow us to clarify some aspects of power-law dis- 
tributions, which have often been overlooked in the past 
literature. It will be shown that the distributions with a 
finite energy cutoff have properties intermediate between 
the microcanonical and the canonical ensembles (which 
are exactly obtained in the two limiting cases q = and 
5=1 respectively) , whereas those with an infinite power 
tail represent a basically new physical situation. The con- 
sideration of the parameter q alone can be misleading. In 
the case of large systems, q is in fact always necessarily 



close to 1, but we will show that this by no means im- 
plies that the character of the distribution is necessarily 
close to Boltzmann-Gibbs. We propose the introduction 
of a new parameter i], expressed as a function of both 
q and the number of degrees of freedom, which provides 
an estimate of the actual deviations from the canonical 
distribution. 

Sections HH - IIVI deal with some basic facts about sta- 
tistical thermodynamics, power-law distributions and 
power-law densities of states respectively. Although some 
of the results we shall obtain in these sections are not 
completely new, we shall present them with an emphasis 
on those aspects which arc at the basis of our particu- 
lar approach to generalized thermodynamics. This will 
lead us in Sec. |V] to our new classification of power-law 
distributions by means of the parameter 77, and in Sec. 
I VII to the main results about the statistical expression of 
thermodynamic temperature and entropy. 

II. PRELIMINARIES 
A. Temperature and thermal equilibrium 

Let us consider a system S composed of two 
weakly coupled subsystems Si and S 2 . We shortly 
denote with z\ the set of 1v\ canonical variables 
5(1,1, . . . , qi. Vl ,PiA, ■ ■ ■ ,Pi,vi of the phase space Ti of the 
subsystem S\, and with Hi[zi) the corresponding hamil- 
tonian function. If the same notation is used for the 
subsystem 1S2, the fact that the coupling between the 
two is weak allows us to express with good approxima- 
tion the hamiltonian H of the complete system S as 
H ~H 1 {z 1 ) + H 2 {z 2 ). 

We define the function Q,i(Ei) as the volume of the 
region of space Ti corresponding to energies of the sub- 
system S\ lower than E\. Introducing the step function 
9 such that 9(x) = 1 for x > 0, 6{x) = for x < 0, we 
can write 

Sl 1 {E l )= [ dz 1 9(E 1 -H 1 (z 1 )) , (1) 

where dz\ stands for a 2i/i-dimensional volume element 
of T\. The derivative U)\{E\) of the above function can 
then be identified with the "density of states" of Si: 

ui(£7i) = £l'i(Ei) = [ dzi 5(Ei - Hi( Zl )) . 

After defining in the same way the functions ft 2 and io 2 
for the subsystem 6>2, we obtain [l?], [HI for the corre- 
sponding functions of the complete system S: 

tl(E) = J dzidz 2 9{E - Hi(zi) - H 2 {z 2 )) 
= J dE\ uji(Ei)fl 2 (E — Ei) 
= J dE 2 Qi(E ~ E 2 )uj 2 (E 2 ) . (2) 
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With each macroscopic state of the system is associ- 
ated a probability density function p on the space T. As 
a consequence of Liouville's theorem [TH, a stationary 
probability density can only be a function of constants of 
motion of the system. It is usual in statistical mechanics 
to assume that p is actually only a function of the en- 
ergy. We will then write p{z\, z 2 ) — 3>(-ffi(zi) + #2(22)), 
where $ is a function of a single real variable, on which 
no particular hypothesis is made for the moment. The 
marginal distribution pi for the subsystem Si, defined 
as pi(zi) = dz2 p(zi, z%), can then be expressed as 
pi(zi) = $i(i/i(zi)), where 

= I dz 2 <S>(E 1 +H 2 (z 2 )) 

= J dE 2 uj 2 {E 2 )<S>{E l +E 2 ) . (3a) 



We have similarly p 2 (z 2 ) = & 2 (H 2 (z 2 )), with 
$ 2 (£ 2 ) = J dE 1 u 1 (E 1 )9(Ei + E 2 ) 



(3b) 



Since u> (E)$(E) represents the probability density 
function for the energy, defining the functional as the 
average of the function Q,(E)/lo{E) we have 



9 = / dE <f>(E)Sl(E) 



(4) 



and similarly for 0i and 02 . Making use of Eqs. @ and 
© we find 

JdE$(E)n{E) = J dEi$i(Ei)fli(Ei) 



dE 2 $ 2 (E 2 )n 2 (E 2 



or equivalently [l7| 



= 0! = 2 . 



(5) 



For the case of a system S described by a microcanon- 
ical ensemble, the quantity was given the name of 
"empirical temperature" in Ref. [2l|, on the grounds of 
the equipartition theorem. For our present purposes, the 
main point is that Eq. ([5]) has here been proved in a com- 
pletely general way, independently of any assumptions 
about the specific form of the probability density or of the 
density-of-states functions. As we explained in the Intro- 
duction, we assume that the parts of an isolated system 
in a stationary state are in thermal (meta)equilibrium 
with each other. We can then affirm that the functional 
assumes the same value for any two systems for which 
this equilibrium condition is fulfilled. Such a result is not 
sufficient to conclude at this point that represents the 
correct generalized definition of thermodynamical tem- 
perature, since the second law has not yet been taken 
into consideration. Besides, there could exist in princi- 
ple other quantities independent of which, at least for 



a certain class of systems, enjoy a property of the same 
form as that expressed by Eq. ([5]) . If we denote the set 
of such hypothetical quantities with r, we can provision- 
ally conclude that a generalized formula for the inverse 
temperature [3 in statistical mechanics, in order to be 
consistent with the zeroth law of thermodynamics, has 
to be expressible as an appropriate function of and r 
of the form 



B. Temperature and thermodynamics 



(6) 



If the hamiltonian H depends, besides on z, also on 
some external parameters aj (such as for instance the vol- 
ume, or the linear dimensions of the system's container), 
then under the infinitesimal variations ddj of these pa- 
rameters the system performs on the surrounding envi- 
ronment the average work 



SW rav = 



dz pdH 



(7) 



with dH = ^2j(dH / da,j)daj . From the first law of ther- 
modynamics we get 



SQre 



dU + SW„ 

d ( I dz pH j — J dz pdH 



dz 5pH 



(8) 



The above equation provides the expression in the lan- 
guage of statistical mechanics for the quantity of heat 
5Q rev which is exchanged by the system in an infinitesi- 
mal reversible transformation. 

According to the so called Clausius theorem, which is 
based on the second law of thermodynamics, 5Q rcv ad- 
mits an integrating factor which can be identified with 
the inverse of the absolute temperature. It can be easily 
shown that a sufficient condition for the validity of an 
equivalent theorem in statistical thermodynamics is the 
existence of a real function F, with a concavity of defi- 
nite sign, such that the probability distribution p can be 
obtained as a stationary point of the functional 



8 = J dzF{p) 



under the constraints 



and 



dz p = 1 



dzpH = U 



(9) 



(10a) 



(10b) 
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Equation (|10a|l simply imposes the correct normalization 
for the probability density, while Eq. (|10b|) fixes the value 
U for the mean energy of the system. Introducing the 
Lagrange multipliers a and /3, such a constrained extreme 
can be located as the zero of a functional derivative in 
the following way: 

= (§-a J dzp-P J dz pH^j 

= F'{p{z))-a-pH{z) . (11) 

The hypothesis we have made on the concavity of F 
guarantees that F' is a monotone function, and there- 
fore admits an inverse function which we call G. From 
the last equation it then follows that we can write p(z) = 
<$>(H(z)) with 

$(E) =G(a + (3E) . (12) 
Equation ([TT]) implies that 

d [s - a J dzp-/3 j dz pH^j = 

for any arbitrary variations 5p(z) at fixed aj. It follows 

dS = ad(^J dzpj+fi J dz SpH = /3SQ ICV , (13) 

where for the last step we have used Eq. ([5]) and the fact 
that d(J r dz p) = on account of Eq. (|10a|) . The above 
equation shows that the state function /3, defined as the 
Lagrange multiplier which solves the system of equations 
(|TU)) - pd|) . is an integrating factor for 5Q 1CV . 

Such an integrating factor is in principle not unique. 
Let us define 

S = kg(S) , (14) 

where k is the Boltzmann constant and g any monotone 
function. We then have, denoting with g' the derivative 
of g, 

dS = kg'(S)dS = k(38Q ICV , (15) 

with 

(3 = 0g'(S) . (16) 

Equation (|15[) shows that also the new state function /3, 
defined by Eq. (fll)|) , is an integrating factor for <5Q rcv . 
The analytical expressions of the functions / and g ap- 
pearing respectively in Eqs. (JBJ) and (fT4"|) are up to this 
point undetermined: they are related to each other by 
the condition that they must provide the same value for 
P. 



C. The Boltzmann— Gibbs distribution 

If we put in Eq. ([9]) F(p) = —p\og(h u p), where h is 
a constant with the dimensions of an action and 1v is 
the dimension of the phase space T, we obtain from Eq. 
(fT2)) $(E) = Z~ 1 e~P E , where the partition function Z = 
h u e 1+a provides the required normalization factor for the 
probability density. From Eq. ([J) we obtain [ItJ 

6 = - / dEe~ 0E n(E) 

= TeU E ^ E)e -"-j <17) 

where, integrating by parts, we have used the fact that 
il(E m i n ) = 0. Hence, according to the considerations of 
the previous subsection we can simply put f r {Q) = 1/0 
in Eq. © and g(S) = S in Eq. (16]), obtaining 

/3=l/6 = /3. (18) 

Furthermore, for an ideal gas composed of a large num- 
ber N of molecules, it is found that (3 = PV/N, where 
P is the pressure and V the volume occupied by the gas. 
Therefore, on the basis of the laws of thermodynamics 
and of the state equation of perfect gases, one can right- 
fully associate with the Boltzmann-Gibbs distribution 

<P(E) = Z- y e-^ E (19) 

the absolute temperature T = \/k(3 and the entropy 

S BG = -k J dzp\og{h»p) . (20) 

It has been assumed until now that the system can 
be correctly described by classical mechanics. However, 
if we identify h with the Planck constant, according to 
Heisenberg uncertainty principle we can associate with 
each element of T of volume h v a quantum state with 
occupancy probability P = h v p. On the basis of this 
correspondence we can rewrite the Boltzmann-Gibbs en- 
tropy (|20[) in the form 

S BG = -kJ2^ogP , (21) 

i 

the sum being extended to a complete set of quantum 
states of the system. 

III. POWER-LAW DISTRIBUTION 
FUNCTIONS 

A. Definition 

Let us fix the arbitrary additive constant of the hamil- 
tonian function in such a way that its minimum value is 
zero. We refer as "power-law distributions" to stationary 
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probability density functions p{z) = $(if(z)) for which 
the function <£>, defined on the interval [0, +00), takes one 
of the two alternative forms 

$(E) = D- 1 (E-E a )- p for < E < +00 (22a) 

with E < 0, or 

HE)= {V-HE -E)- P for 0<E<E 
v ' |0 for £>£ . y ' 

In the above equations p is a real parameter, and D a 
normalization coefficient such that 



+ OC 



dEuj(E)<&{E) = 1 



(23) 



In Eq. (|22b[) the usual Tsallis cut-off prescription @, [HJ 
has been adopted. It has become customary in the litera- 
ture to equivalently express these distributions by means 
of the ^-exponential function [2j , but this would not here 
be helpful in view of the formal manipulations we shall 
need to perform. 

The requirement that the integral on the left-hand side 
of the last equation be finite imposes limitations on the 
possible values of p. In the case of Eq. (|22a|) , if we assume 
that the behavior of the density-of-states function in the 
limit of high energies is expressible as 



u(E) oc E s 



for E 



(24) 



with s > 0, then a necessary condition is p > s. Fur- 
thermore, if the first law of thermodynamics has to hold, 
it must be possible to define the average energy of the 
system as 



U = 



dEoj(E)<f>(E)E 



(25) 



Therefore, in order for U to be finite, in the case of Eq. 
(|22a[) it is actually necessary that p > s + 1. In the 
case instead of Eq. (|22b[) . the requirement that both the 
integrals in Eqs. (j2"3")) and (|2"5} be convergent in a left 
neighborhood of Eq is equivalent to the condition p < 1. 

In the rest of the present paper we shall often write 
the power-law distributions in the general form 



=D- 1 \E-E \ 



(26) 



where it is understood that E can only take values such 
that E-E > (resp. E — E o <0) whenever p > 1 (resp. 
p < 1). The qualitative behavior of the function <!> in the 
four cases p>l,0<p<l,— I < p < and p < — 1 is 
displayed for the sake of clarity in Fig. 1. When Eq. (|24p 
holds, for the reasons we have explained p can never lie 
within the interval [1, 1 + s]. 



B. g-Distributions and thermodynamics 

It is easy to see that a power-law distribution can 
be obtained via a maximization procedure like that il- 
lustrated in Sec. Ill Bl In fact, if in Eq. © we take 
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FIG. 1: Qualitative behavior of the function $ for different 
values of p. 



F(p) = h u ( q ^p q , where q is a real number, then Eq. 
(fT2l) provides 



$(E) = hr 



a + /3E 



1/(9-1) 



which coincides with Eq. (|26j) when putting 

1 



P = 



(27) 



(28) 



E = (29) 
D = h^p/ql 1 /^ . (30) 

Looking at Eq. (|27p one also sees that $ will be of the 
form (|22ap when 0/q > 0, and we have shown that in 
such a case the existence of a finite mean energy requires 
p > 1 as a necessary condition. Viceversa, $ will be of 
the form (|22b[) when 0/q < 0, and in that case we must 
then have p < 1 . 

We can rewrite Eq. ([9]) in the form 



s q = h~ v J dz (h u P y , 



(31) 



which in terms of the discrete probabilities of the quan- 
tum states becomes 



(32) 



According to Eq. ([14]) the entropy will be of the general 
form 



Sq — kgq(Sq) 



(33) 



where g q is an appropriate monotone function which may 
contain q as a parameter. Tsallis originally derived his 
power-law distribution from the function [7J] 



2-1 



(34) 
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On the other hand, the expression 

}ogS q 



S q - fc-j 



9 



(35) 



had already been considered by Renyi in the context of 
information theory 19]. Other variants have also occa- 
sionally received some attention in the literature: we may 
here mention Tsallis normalized entropy [23| 



S» = k 1 Sq . 
and Tsallis "escort" entropy [lij 



< /. '/ 1 ^ :7 



1-9 



(36) 



(37) 



Note that Pf ~ P 4 + (g - l)piogP for <? ->• 1, whence 
5 g ~ 1 + (q — 1) J2 t Pi log P^. It is then easy to see that 



all the four functions given in Eqs. (|3"4"]) - (f3"T]) tend to the 
Boltzmann-Gibbs entropy (j2Tj) for q — > 1. 

From the considerations made in Sec. Ill Bl it is clear 
that all these functions (and other possible ones) deserve 
equally well the name of entropy as far as the second 
law of thermodynamics is concerned, since for each of 
them one can define a state function such that Eq. 
([TS"]) is satisfied. The criterium for the correct choice is 
provided in fact by the zeroth law. In the following of 
the present paper, after a closer analysis of the condition 
of thermal equilibrium between power-law distributions, 
we shall aim at determining the analytic form of g q for 
which can be expressed in the form ©, with 9 given 
by Eq. flU). 



C. A generalized partition function 

When a system is described by a power-law distribu- 
tion of the form ([26"]) , some of the quantities previously 
introduced can be conveniently expressed with the aid of 
the auxiliary function 



C(E ) 



dEuj{E)\E - E 



i-p 



1 1 - p\ J dE Q(E) | E - Eq \- p 



(38) 



which, in the same way as an analogous one introduced in 
Ref. [15[ , plays the role of a generalized partition function 
for this class of systems. According to the considerations 
made in Sec. IIII Al when p > 1 the function £ is defined 
for Eq < and the integration domain / in Eq. ([38]) is 
the interval [0, +oo), whereas when p < 1 we have Eq > 
and I = [0,Eo). 

From Eq. we get 



D(E ) = / dEuj{E)\E - Eq 



C(Eq) 



(39) 



We can then express the function O defined by Eq. (|4|) 



as 



9 = Dk)L dEnmE - Eol ~ p 

C(E ) 



\l-p\D(E ) 



dE 



logC(P ; 



(40) 



We can also evaluate the mean energy U using the 
relation 

\U-E \ = J ' dEw{E)${E)\E -E Q \ 



CCgo) 

D{Eq) 



= |i-p|e 



(41) 



Since it is obviously U > Eq (resp. U < Eq) for p > 1 
(resp. p < 1), Eq. (|4"Tj) is equivalent to 



u = E + (p-i)e . 

From Eq. (f3Tj) we finally obtain 



(42) 



~ s ^=wShL dEu)mE - Eorvq 

whence, using Eqs. H2SJ), (001) and (gTJ), 

S q = h- v l*>[D(E Q )]Vv- l (;{EQ) 
C(Eo) _ 



D(E ) 1 - q 



. 



(43) 



The expressions obtained in this subsection will be use- 
ful in the following of the paper, when we will have at 
our disposal an explicit expression for the function £. 



IV. POWER-LAW DENSITY-OF-STATES 
FUNCTIONS 

Let us suppose that the functions f2 of the two weakly 
interacting systems Si and S2, already considered in Sec. 
Ill Al are given respectively by Qi(Ei) = BiE^ 1 and 
D, 2 (E 2 ) = B 2 E 2 2 , with s ± > 0, s 2 > 0. According to 
Eq. @ the function f2 for the system S, composed of 1S1 
and 1S2, is then given by 

n(E) = PiP 2 si / dEi E{ 1 ^ 1 {E - El) 32 
Jo 

= B 1 B 2 s 1 E Sl+S2 [ dxx^-^l-xY 2 . 
Jo 

Since the integral in the last expression has the value 
r(si)r(s 2 + 1)/T(«i + s 2 + 1) SIS, we obtain 



Q(E) = BE S 



(44) 
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with 



and 



s = Si + s 2 



(45) 



B — BiB 2 - 



r ( 8l + i)r(s 2 + i) 
r(si + s 2 + 1) 



Note that the last equation can also be put in the expres- 
sive form 

\og[BT{s + 1)] = log[Sir(si + 1)] + log[S 2 r(s 2 + 1)] . 

(46) 

From Eq. flUJ) it also follows that u(E) = BsE s ~ l . 

Let us further assume that the system S is described 
by a power-law probability distribution of the form (|26[) . 
The function £ defined by Eq. ([55)1 can be explicitly cal- 
culated as 

C(JSb) = Bs J dEE 3 - 1 \E-E \ 1 - p = A\E \ s - p+1 , (47) 
where [HJ for p < 1 

io r( s - P + 2) 

(48a) 

whereas for p > s + 1 (recall that values of p between 1 
and 1 + s are forbidden) 



A = Bt 



dx 



o {x + iy- 1 



= B 



r(s + i)r(p-*-i) 



r(p-i) 



By using Eq. ([55]) we then obtain 



p-l 



(48b) 



(49) 



We can also calculate the marginal distribution for the 
subsystem S\ using Eq. (f3"aj) . This gives 

= ^jy-f j dE 2 E s ^- 1 \E 1 +E 2 -E \- p 

where h — [0, +oo) for p > 1 and /i = [0, E — E\) for 
p < 1. It is then immediate to see that 



1> 1 (E 1 ) = D^ 1 \E 1 -Eo\- pi , 



(50a) 



where pi = p — s 2 , while the coefficient D\ can be ex- 
pressed as a function of Eq, B\, S\, and pi according to 
the same formulae (|49p and ([4"B")) which give the depen- 
dence of D on _Bo, -B, s, and p. In the same way, from 
Eq. J3b} one finds 



<f> 2 {E 2 ) = D^ 1 ^ - E \- p2 



(50b) 



with p 2 = p — s±. When p < 1 it is also obvious that 
^ X (E) = $ 2 (J5) = for B > S . 

Equations ([5^|) show the remarkable fact that, when 
the densities of states of the systems under consideration 



exhibit a power-law dependence on the energy, power-law 
probability distribution functions are compatible with 
the condition of thermal equilibrium which was described 
in Sec. Ill Al We have in fact that the two weakly in- 
teracting systems iSi and S 2 , and the composed system 
S, are all simultaneously described by stationary power- 
law distributions of the form . We note also that all 
these distributions have in common the same value of Eq , 
whilst the relation among p\, p 2 , and p can be written as 

r = n = r 2 , (51) 

where we have introduced the new parameter 



r = p — s 



(52) 



This means that Eq and r are intensive quantities which, 
like 0, assume the same value for two systems at mu- 
tual equilibrium, independently of the respective size and 
physical nature. From Eqs. (|4U|) and ||4"T|) we obtain 
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En 



1 



(53) 



which implies that only two of these three quantities are 
actually independent. It follows that we can identify r 
with the quantity that was introduced as a parameter of 
the function f r in Eq. (J6j> . On the other hand, Eqs. |45|) 
and (|4*6")1 show that the quantities s and log [Br (s + 1)] 
are extensive, since their values for the complete system 
S are the sum of the corresponding values for the two 
subsystems Si and S 2 . The same conclusions can obvi- 
ously be extended to an arbitrary number of subsystems. 

For the reasons that we have illustrated, systems pos- 
sessing both a power-law density-of-states and a power- 
law probability distribution function (at least within an 
acceptable approximation) are likely to be the most sig- 
nificant ones with respect to their thermodynamic prop- 
erties. In the following sections we shall refer to them 
shortly as "power-law systems" . 



A CLASSIFICATION OF POWER-LAW 
SYSTEMS 



Let us consider a power-law system characterized by 
Eqs. gl} and Then from Eqs. @2J| and {53} we 

obtain for the mean energy 



U 



s — p + 1 



E = sQ . 



(54) 



In order to make a comparison with the Boltzmann- 
Gibbs distribution, let us rewrite the function $ in the 
form 

^(E)=D- 1 \E \-P e xp[-plog(l-E/E )] . (55) 

For E < \E \ we have log(l - E/E ) = -E/E + 
0(E 2 1 Eq), and the term quadratic in the energy gives 
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a negligible contribution to the exponential function in 
Eq. (|55p whenever 



(56) 



It is now clear that, in order to evaluate the overall be- 
havior of the statistical ensemble under consideration, 
only energies of the same order of magnitude as U need 
be considered, since the probability that the system may 
take energy values much higher than U is in general van- 
ishingly small. If we put E = xU = E xs/(s — p + 1), 
where a; is a dimensionless variable, we see that the in- 
equality E <C \Eq\ will be satisfied for all x up to order 
unity provided that \p— s — 1| 3> s, whereas the condition 
([56]) becomes equivalent to (p — s — 1) 2 /H 3> s 2 . It is 
then easy to see that these two conditions will both hold 
only when \p\ ;§> s 2 . Recalling Eq. Q28p . this is equivalent 
to 

|«-i| = n«-?- ( 57 ) 

\P\ s 

In this limit we can rewrite Eq. (|55[) as 

~ D-^Eol-'expipE/Eo) , (58) 

which corresponds to a canonical distribution with pa- 
rameter f3 = -~p/Eq. We thus find again the well-known 
result that a q-distribution tends to the Boltzmann- 
Gibbs distribution in the limit q — ► 1. In addition, the 
above analysis gives the size of the neighborhood of 1 in 
which the parameter q has to fall in order for the two dis- 
tributions to be practically equivalent over the relevant 
range of energies. This size is of order 1/s 2 , and thus 
decreases as the inverse square of the dimension of the 
system (recall for instance that, for an ideal monoatomic 
gas consisting of N particles, one has s — 3 AT/2). 

We already know from Sec. IIII Al that for < p < 1 
one has <S>(E) = for E > E and 

lim $(E) = +oo . 

E^E- 

Recalling Eqs. (|48a|) and 1(15]) we can write for E < E Q 



$(£) 



T(s-p+l) 
BT(s + l)T(l -p) 



Et s {E Q -E)- 



whence 



HE) 



1-p 



BsE s -\E -E) 



for p 



(59) 



Therefore <fr(E) vanishes in the limit p — > 1 _ for any 
E =/= Eq, and the probability density on the phase space 
becomes fully concentrated on the surface of equation 
H(z) — Eq. This means that the power-law distribution 
tends to the microcanonical ensemble in the limit p — > 
1 _ , which in turn corresponds to q — > _ . If we consider 
again a system S with a subsystem Si, and we suppose 



that S has a microcanonical distribution, so that p ~ 1, 
then using Eq. (p)Tj) we havepi = p— s + si < 1, since s > 
S\. This is in agreement with the observation, already 
made some years ago [ll[ , that a power-law distribution 
with finite energy cutoff may arise as the marginal of a 
microcanonical distribution for a finite system. 

The region s + 1 < p < s 2 is the one in which power-law 
distributions display features which are most remarkably 
different from other known ensembles. Here the tail of 
the distribution at high energies becomes relevant, and 
represents a significant departure from the exponential 
decrease of the canonical ensemble. The consequences 
of this phenomenon become extreme in the limit p —> 
(s+l) + , which corresponds to q — > [s/ (s+l)] + . Note that 
for macroscopic systems, having very large s, this limit 
value for q still looks extremely close to 1, i.e., to the value 
which characterizes the Boltzmann-Gibbs distribution. 
Therefore q does not appear to be the right parameter 
to distinguish among the different possible behaviors of 
a power-law system. To this purpose it is instead useful 
to introduce the parameter 



P-1 
p — s — 1 



= 1 



T - 1 



(60) 



Recalling that p can vary in the set (— oo, l)U(s+l, +oo), 
we see that rj takes values in the full interval (0,+oo). 
The canonical ensemble, for which q = 1, p = ±oo, cor- 
responds to rj = 1; the microcanonical limit q —> CP, 
p — > 1~ corresponds to rj — > + , and the "long power- 
tail" limit p — > (s + 1) + , q — > [s/(s + 1)] + corresponds to 
r/ ^> +oo. For < r] < 1 the power-law distribution is in 
some sense intermediate between the microcanonical and 
the canonical ensembles. According to the previous dis- 
cussion, the behavior of the system becomes essentially 
canonical for \q — 1| <C 1/s. A scheme of the behavior of 
the parameters p, q and rj is reported in Table |TJ 

Let us write down the relation between the parameters 
r\ and r\\ of two systems S and Si at equilibrium with each 
other. Since r — 1 = s/(rj — 1), from Eq. (|5ip we obtain 



si 



n - 1 



1 + Si- 



(61) 



If we now consider a large power-law system S composed 
of many weakly interacting microscopical particles, it is 
interesting to take as Si the subsystem constituted by 
a single particle, for which the parameter si is typically 
of order unity (for an ideal monoatomic gas one has for 
instance s = 3N/2, si = 3/2). According to the previous 
analysis, the one-particle distribution will be appreciably 
different from Boltzmann-Gibbs when \r)\ — 1| > 1/si = 
O(l), which on account of Eq. (|6"Tj) is equivalent to \rj — 
l| s i/ s > 0(1). For s 3> si this can clearly occur only 
when rj > s, which corresponds to a system S far in the 
long power-tail regime. 
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TABLE I: Values of the parameters p, q and rj for the different possible types of power-law systems. Symbols f and \ indicate 
increase or decrease respectively. 





microcanonical 






canonical 


long power-tail 


p 


1 


\ 


x 


s, -oo;+oo N 


V s + 1 


q 





\ 


— oo; +oo x 


1 


%, s/(s + l) 


V 





/ 




1 


^ +00 



VI. THERMODYNAMIC FUNCTIONS OF 
POWER-LAW SYSTEMS 

A. Temperature and entropy 

Let us now consider the crucial problem of determining 
the thermodynamic temperature and entropy of a power- 
law system. According to Eqs. ^ and (fTo| we must have 



= f r (Q) = pg q (S q ) 



(62) 



Taking into account Eq. (|43|) , the above equation can be 
rewritten in the form 



9/ r (6) = (1 - q)S q g'(S q ) . 



(63) 



When we introduced the function g q (S q ) in Eq. (f3"3"]l . we 
made the assumption that it may contain the only param- 
eter q (see the last part of Appendix [X] for some related 
considerations). We then observe that the two variables 
r and 8, that appear on the left hand side of Eq. (p)3"l). are 
independent of the two variables q and S q that appear 
on the right hand side. It can in fact be easily checked 
(see Appendix [5} that all these four quantities can be 
expressed as independent functions of the four free pa- 
rameters p, Eq, s, and B that characterize the system 
according to Eqs. (|2"6"|) and (|4"4")l . It is then obvious that 
the two members of Eq. (|63p can be identically equal to 
each other only if they are both equal to a common con- 
stant value c. If we first consider the second member, we 
can thus write 



9a(S q ) = 



(i - i)s q 



(64) 



Integrating this equation and recalling Eq. ([33]) . we ob- 
tain that the entropy S q of the power-law system is given 
by 



kg q (S q 



ck 



\ogS q 



1 



1(9) . 



(65) 



where I(q) is an a priori arbitrary function of q. The 
requirement that the Boltzmann-Gibbs expression (|2ip 
has to be recovered in the limit <7 — > 1 imposes the con- 
ditions c = 1 and 1(1) = 0. If we further require that, 
in accordance with the third law of thermodynamics, the 
entropy vanishes at zero temperature for all systems with 



a nondegenerate ground state, we must set I(q) = for 
any q. In this way we finally obtain 



So 



'l-q 



(66) 



which coincides with the expression of Renyi entropy Sjf 
given by Eq. l[3"5"]). 

Setting also the left-hand side of Eq. (f6"3"|) equal to 1 
provides the expression for the inverse absolute temper- 
ature: 



1 

kf 



/3 = /r(6) 



1 

e 



(67) 



We therefore find that the relation kT = 9, which had 
already been proved for the Boltzmann-Gibbs statistics, 
is valid for all power-law systems. From Eqs. (|53[) and 
([54")) we also find for the mean energy U of the system 
the value 



skT 



(68) 



As in standard thermodynamics, the energy is thus di- 
rectly proportional to the absolute temperature T, and 
the specific heat at constant volume is given by 



Cv = ^ = sk 



(69) 



As a verification of the results obtained in this section, 
we examine in Appendix [B] the consequences of adopting, 
instead of Eq. ([35]) . the three different entropies men- 
tioned at the end of Sec. IIIIBI 



B. Properties of Renyi entropy 

On account of the results of the previous subsection, 
from now on we shall directly write in place of S q . 
From Eqs. (J43J) , flU, (gSJ) and j53]) we obtain 



k~ S. 



1 qR 



-vlogh + log A + (p — 1) log 
- slogflp - s - l\kT) . 



p — s — 1 



(70) 



Having in mind that the number of degrees of freedom 
v and the value of the external quantities such as the 
volume are directly related to the parameters s and B of 
the density of states (|4"4"|) , and recalling Eq. (|2"8|) , we want 
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here to analyze the behavior of the entropy considered as 
a function of the independent variables q, a, B, and T. 

Let us first of all compare Eq. ([70| with the corre- 
sponding expression for the Boltzmann-Gibbs entropy. 
By substituting for p = &(E) the expression (p~9|) into 
Eq. ([20]) . and observing that, on account of Eq. (1441) . 



dze~' 3H = / dE sBE s ~ 1 e~P E 
r Jo 



Z 



one obtains 



k^SBcis-B.T) = -v log ft+log[Br(s+l)]+s log kT+s . 

(71) 

We can write in general 

S*(s,B,T) = Sbg(s,B,T) + AS q (s) , (72) 
where, using for A the expressions (|48p . we have 

fc- x AS,(a) = log r , r( l7 P) S +Plog(l~p) 
T(s + l—p) 

+ (s-p) log(a + l-p)-s (73a) 

for p < 1 , and 

k^AS^s) = log r ^~ s) + plog(p-l) 

- (p - a) log(p - a - 1) - a (73b) 

for p > a + 1 . By calculating the limit of the expression 
(|73a[) for p — > — oo and the limit of (|73b[) for p — » +oo 
one can easily check that lim,^i AS q (s) = 0. We can 
therefore conclude that, as expected, 

lim S?(s,B,T)=S BG (s,B,T) . 

Another interesting case to consider is the limit q — > 
0~, or p — > 1~, which according to our discussion of Sec. 
fVl corresponds to a microcanonical distribution. Taking 
into account Eq. we obtain from Eqs. (fTTj) and (|73a|) 

lim S*(s,B,T) = k(-p\ogh + logB + slogU) 



k log 



ft (17) 



(74) 



which is indeed a well-known expression for the entropy 
in the microcanonical ensemble [IT], [HI • 

Since the expression (|7ip is manifestly extensive, in 
order to investigate the thermodynamic limit it is only 
necessary to study the behavior of AS q for a — > +oo. If 
we treat 5 (and therefore also p) as a constant, recalling 
that p can never lie in the interval [1, s + 1] we see that 
such a limit can only be taken when p < 1. Therefore, 
applying Stirling's formula to the expression (|73a[l we 
obtain 



From a physical point of view it might seem more ap- 
propriate to consider a thermodynamic limit in which 
the intensive quantity r is kept constant instead of q. 
We have in this case p = a + r — * +00, which implies 
p > s + 1. Then from Eq. (|73b|) we obtain 



AS, (a) 



ifclogs 



for 



(76) 



with r — constant, q = 1 — l/(s + r) — > 1 . 

Let us finally consider the thermodynamic limit in 
which we keep constant the parameter r\ introduced in 
Sec. El We have in this case to substitute p = rjs/(r] — 
1) + 1 into Eq. (j73a|) when < rj < 1, or into Eq. (|73b|) 
when rj > 1, and then to take the limit for s — > +00. We 
obtain in this way 



lim AS q (s) = 



->+oo 



§ log 77 
-ifclog?? 



for 
for 



< i] < 1 
77 > 1 



(77) 



AS, (a) 



-fclogs for s — > +00 . 



(75) 



with 7/ = constant, g = a/(a + 1 — 77 ) — > 1. 

We see that in the limits considered AS, either in- 
creases logarithmically with the size of the system, as in 
Eqs. (fT5"]) - (fT6"]) , or tends to a constant value as in Eq. 
(|77|) . It follows that in all cases the extensive part Sbg 
prevails in Eq. ([72"]) . and the entropy S^(a, B,T) essen- 
tially reduces to Sbg(s, B,T) as s -» +00. We note 
also that q does not appear in the expression (f6"5)) of the 
mean energy. These results can be seen as a general- 
ization of the well-known theorem about the equivalence 
of the different ensembles in the thermodynamic limit. 
We have shown in fact that this equivalence does not 
only hold between the canonical and the microcanonical 
ensembles, but also among the power-law distributions 
corresponding to all possible values of rj. This fact is 
all the more remarkable if one recalls that, as we have 
illustrated in detail in Sec.|Vl the microscopical distribu- 
tion function can instead be significantly different from 
Boltzmann-Gibbs also for very large systems. 

We would like finally to mention that the additivity 
of Renyi entropy for mutually independent systems was 
already well known. Our analysis has shown on the other 
hand that it is also additive for large power-law systems 
which, instead of being independent, are at mutual ther- 
mal equilibrium in the sense that we have specified in 
Sec. Ill Al One can observe nevertheless that the Renyi 
entropy (|72[) includes, at variance with the Boltzmann- 
Gibbs one, a nonextensive contribution AS, which may 
be nonnegligible for systems of small size. 



C. The ideal gas 

Let us consider as a significant example the case of an 
ideal monoatomic gas composed of N particles of mass 
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to in a volume V. We have 



v 
s 

B 



3N 
3N 

{2irmf N ' 2 V N 
T(3iV/2 + 1)N\ 



where the factor N\ in the denominator of the last expres- 
sion accounts for the indistinguishability of the atoms. 
From Eqs. (|68|) and (jTTj) we then obtain for large N 



U 



-NkT 



Si 



2 

Nk 



3 2nmkT 
2 bg ^^ 



log 



V_ 

N 



OQagN) 



(78) 



(79) 



We find in this way two formulae which are already famil- 
iar from standard statistical thermodynamics. It is then 
easy to see that ideal gases governed by power-law dis- 
tributions obey the same state equation of Boltzmann- 
Gibbs ideal gases. Consider in fact an infinitesimal 
isothermal transformation, for which dT = 0. Then 
from Eq. ([75)1 we get dU = 0, and the infinitesimal work 
SW = PdV performed by the system, where P is the 
pressure of the gas, can be obtained according to the 
first and second laws of thermodynamics as 



PdV = SQ = TdSj; = NkT 



dV 
~V 



(80) 



where the last equality follows from Eq. (|75j) . Compar- 
ing the first and last members of the above equation we 
finally get 



the main new ingredient, with respect to former inves- 
tigations, is the observation that Tsallis statistics is ca- 
pable of describing stationary states, in which power-law 
systems combine together in such a way to give rise to 
a larger composed power-law system. The consideration 
of such a circumstance allows one to extend to a deeper 
level the analogy with the standard thermodynamics of 
weakly interacting systems [171 Il8j . which has already 
proved to be an illuminatin g gu ideline for the analysis of 
Tsallis thermostatistics [OJll, El, E3, El, El] • 

Of course power-law distributions can in principle ap- 
ply also to systems for which the hypotheses of weak 
coupling or of power-law density of states are in general 
not satisfied. An important case of this type is likely to 
be represented by systems of particles with long-range in- 
teractions, which are indeed considered among the most 
interesting candidates for the applicability of nonexten- 
sive thermodynamics. In such cases the conclusions of 
our analysis may not necessarily hold, and it is possible 
that other theoretical investigations are to be carried out 
on the basis of a completely different approach. In par- 
ticular, other descriptions might become appropriate, in 
which no reference at all is made to the usual formalism 
of equilibrium thermodynamics. For instance, a justifi- 
cation for Tsallis entropy has recently been proposed on 
the basis of purely dynamical assumptions about the be- 
havior of a system far from equilibrium [271 ] . It has also 
been pointed out that deviations from the exact Tsal- 
lis distribution law can be present in real situations [9(. 
It is reasonable to expect that such deviations may be 
related to the specific form of the density-of-states func- 
tion, so that a generalization of Tsallis statistics might be 
required in order to treat the cases in which the density 
of states significantly departs from a power law. 



PV = NkT 



(81) 
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In particular, one can conclude that also for power-law 
systems the ideal gas temperature defined according to 
Eq. (|8ip is the same as the thermodynamic temperature 
based on the second law. 



VII. CONCLUSION 

We have developed a new generalized thermodynam- 
ics for systems governed by Tsallis distributions, in a way 
which also allows for a rigorous statistical interpretation 
of the zeroth law. Our approach thus provides a new 
point of view on an issue which has been the object of 
sharp debate in recent literature (2r| . The main outcome 
of our analysis is that the thermodynamic entropy associ- 
ated with these distributions is expressible as a function 
of the probabilities according to the formula first intro- 
duced by Renyi. 

Our results have been obtained under a minimum set 
of simple and apparently plausible hypotheses. Probably 
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APPENDIX A: ON THE INDEPENDENT 
PARAMETERS OF A POWER-LAW SYSTEM 

The deduction of Eq. from Eq. (p)3"|) is based on 
the mutual independence of the four variables r, O, q and 
S q . In order to elucidate this point, let us start from the 
four main free parameters which can be put at the basis 
of the description of a power-law system. These are s, B, 
p and Eq. The first two of them appear in the expression 
(|44[) which gives the primitive function of the density of 
states, and are therefore directly related to the form of 
the hamiltonian function. Their physical meaning is il- 
lustrated by the example of the ideal gas which is treated 
in Sec. IVI CI In particular, 2s expresses the number of 
degrees of freedom contributing to the hamiltonian, and 
is therefore related to the size (number of particles) of 
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the system, while B contains the dependence on the ex- 
ternal parameters which, like the volume V, are varied 
when the system exchanges mechanical work with the 
surrounding environment. On the other hand, p and Eq 
are the two free parameters in the power-law probability 
distribution (126|) . D being a normalization constant. Ac- 
cording to Eqs. (f28|) and ([54]) . p is simply a function of 
q, while Eq fixes the mean energy, and thus the temper- 
ature of the distribution. For the case of the ideal gas, 
for instance, the mutual independence of s, B, Eq and p 
amounts to the mutual independence of N, V, T and q. 

Assuming for simplicity the relation v — 2s as for the 
ideal gas (the particular relation between these two pa- 
rameters is inessential with respect to the present con- 
siderations), we can express q, r, and S q as functions 
of the independent quantities s, B, Eq and p, according 
to the relations 



1 



e 



s„ = 



i 

p 

p — s 

Eq 
1 + s - p 

s + 1 - p 



BG(s,p)T(s + l) 



\Eq 



(Al) 
(A2) 
(A3) 

(A4) 



Equations (|A1[) - (|A3|) are identical respectively to Eqs. 
(|25]) . ([52]) and ([53]). Equation l[A"4]) . in which we have 
introduced the function 



G(s,p) 



T(l-p)/T(s-p+l) for p<l 



r(p - s)/T(p) 



for p > s + 1 



is derived from Eqs. @3J and (j4T]) - (|4l)j) . 

It easy to see that the four functions defined by Eqs. 
(|Aip ~ (|A4p are mutually independent. They can in fact 
be explicitly inverted, so as to express s, p, Eq and B as 
functions of q, r, and S q . We obtain 



1 



1 



P 

Eq 

B 



(l-r)6 

K(g,r) 
r(l/(l - q) - r + 1) 

(g-l)(l- r )~ V /(1_9) 



(A5) 

(A6) 
(A7) 



81-r 



l/(l- 9 )-r 



(A8) 



where 



(r(l - r) 
— ^- for q < 0, g > 1 
r(g/(g-i)) 
r 1/ (1-9 



I r(i/(i- g )) 
I r(r) 



for 



8 + 1 



< q < 1 



The argument we have used after Eq. (|63[) , in order to 
justify the introduction of the constant c, is also based on 
the assumption that the function g q (S q ) does not contain 
any parameter other than q. Note that all the functions 
appearing on the right-hand sides of Eqs. ([34]) ~ ([37]) do 
agree with this general assumption, which is dictated by 
obvious reasons of simplicity and economy. We expect 
in fact that such a fundamental physical quantity as the 
entropy should have the simplest analytical expression 
which is compatible with the laws of thermodynamics, 
without the inclusion of unnecessary parameters. The 
argument about the constancy of c is also at the basis 
of the remarkable equality, which we have found in Sec. 
IVI A| between the "empirical temperature" 8 and the 
thermodynamic temperature kT. It is however interest- 
ing to examine the consequences of restraining from the 
above assumption, and considering in Eq. (|33[) a func- 
tion g q ^ r (S q ) which may also depends on the parameter 
r. This would imply the replacement of the constant c 
on the right-hand side of Eq. (|65[) with a function c(r) , 
which would then remain as a multiplicative coefficient 
also in the final expressions of the entropy and of the 
inverse temperature: 



S q = fcc(r) 



1-fl 



T 



e 



fcc(r) 



Taking into account Eqs. (|A2[) and (|A6|) . one should then 
be forced to require that limr^ioo c(r) — 1, in order for 
the Boltzmann-Gibbs entropy and temperature to be re- 
covered in the limit q — > 1. Furthermore, we have shown 
m Sec. E] that a power-law distribution tends to the mi- 
crocanonical ensemble for q — > 0. If one then requires 
that the correct microcanonical entropy ([74]) and the em- 
pirical temperature kT = = U / s [21( be obtained for 
any s > in the limit g — > 0, it would also be necessary 
to impose c(r) — 1 for all r such that — oo < r < 1. These 
can be considered as additional arguments in favor of the 
natural choice of taking c = 1 for all r. 



APPENDIX B: COMPARISON AMONG 
ALTERNATIVE ENTROPIES 



In Sec. IVI Al we have proved that the condition of com- 
patibility with the zeroth law of thermodynamics uni- 
vocally determines the expression of the entropy for a 
power-law system. As a verification of this result, it is 
interesting to check directly how this condition of com- 
patibility is violated by the three particular functions 
which have been mentioned at the end of Sec. IIII Bl as 
well-known possible alternatives to Renyi entropy. 

It follows from the first and second laws of thermody- 
namics that the inverse of the thermodynamic tempera- 
ture (3 must be related to the entropy S by the equation 



k(i 



dS_ 

dU 



(Bl) 
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where the partial derivative of the entropy is evaluated 
at constant external parameters (i.e., those parameters, 
such as the volume of the system, whose variation is re- 
lated to the production of mechanical work). When the 
entropy has the form (|33[) we have then 

P = 9' q (S q )^-. (B2) 
From Eqs. (|i3)l. (|47|) . (l49j) and (54} we obtain 



h~ v A 



p — 1 \ p 1 f \p — s — 1| 



p — s — 1 



C7 



i/p 



According to Eq. the quantity ^4 depends on the 

external parameters (through B) but not on the internal 
energy U . We therefore obtain 



dS n 



dU P U pQ ' 

Substituting for g' q in Eq. (|B2[) the derivatives of the ex- 
pressions which appear on the right-hand sides of Eqs. 
(|34]) - l|3"7[) . we then find for the function (3 in the four 
respective cases 



We thus see that only (3 R , derived from Renyi entropy, is 
equal to the inverse of the parameter which was intro- 
duced in Sec. Ill Al on the basis of the zeroth law of ther- 
modynamics. On the contrary, (3 T , /3 N and /3 s , which are 
respectively associated with the entropies (p4|) . (j36|) and 
(|37| . depend also on S q , which in general takes differ- 
ent values for systems at thermal equilibrium with each 
other. It is therefore confirmed, by the analysis of these 
examples, that only Renyi entropy leads to an expression 
for the thermodynamical temperature which is compat- 
ible with the zeroth law. This result obviously does not 
exclude that Tsallis or other entropies can nevertheless 
be conveniently employed whenever some of the starting 
hypotheses of the present work (such as weak coupling or 
power-law density of states) are manifestly violated by 
the system under investigation. In such circumstances 
a proper modification of the standard laws of thermo- 
dynamics will also be necessary, as it has already been 
recognized in the literature. 



/3 T 
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